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Abstract. Our understanding of Pulsar Wind Nebulae (PWNe), has greatly improved in the last years thanks to unprecedented 
high resolution images taken from the HUBBLE, CHANDRA and XMM satellites. The discovery of complex but similar inner 
features, with the presence of unexpected axisymmetric rings and jets, has prompted a new investigation into the dynamics 
of the interaction of the pulsar winds with the surrounding SNR, which, thanks to the improvement in the computational 
resources, has let to a better understanding of the properties of these objects. On the other hand the discovery of non- 
thermal emission from bow shock PWNe, and of systems with a complex interaction between pulsar and SNR, has led to 
the development of more reliable evolutionary models. I will review the standard theory of PWNe, their evolution, and the 
current status in the modeling of their emission properties, in particular I will show that our evolutionary models are able to 
describe the observations, and that the X-ray emission can now be reproduced with sufficient accuracy, to the point that we 
can use these nebulae to investigate fundamental issues as the properties of relativistic outflows and particle acceleration. 
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INTRODUCTION 

Pulsar Wind Nebulae (PWNe) are bubbles of relativis- 
tic particles and magnetic field created when the ultra- 
relativistic wind from a pulsar interacts with the ambi- 
ent medium, either SNR or ISM. The best example of a 
PWN, often considered the prototype of this entire class 
of objects, is the Crab Nebula. The first theoretical model 
of PWNe was presented by Rees & Gunn [74], developed 
in more details by Kennel & Coroniti [52, 53] (KC84 
hereafter), and is based on a relativistic MHD descrip- 
tion. The ultra-relativistic pulsar wind is confined inside 
the SNR, and slowed down to non relativistic speeds in a 
strong termination shock (TS). At the shock the toroidal 
magnetic field of the wind is compressed, the plasma is 
heated and particles are accelerated to high energies. A 
bubble of high energy particles and magnetic field is pro- 
duced where the post-shock flow expands at a non rela- 
tivistic speeds toward the edge of the nebula. 

Despite its simplicity the MHD model can explain 
many of the observed properties of PWNe. Acceleration 
at the TS accounts for the continuous, non-thermal, very 
broad-band spectrum, extending from Radio to X-rays, 
with spectral index in the range 0-1 .2, steepening with in- 
creasing frequencies, and modeled as synchrotron emis- 
sion |[83,|5l|92,[93l[6^. The under-luminous region, cen- 
tered on the location of the pulsar, is interpreted as the 
ultra-relativistic unshocked wind. Polarization measures 
10, S S S S show that emission is highly polar- 
ized and the nebular magnetic field is mostly toroidal, 
as one would expect from the compression of the pulsar 
wind. This model also predicts that PWNe should appear 
bigger at smaller frequencies: high energy X-rays emit- 



ting particles have a short lifetime for synchrotron losses, 
and they are are present only in the vicinity of the TS; in 
contrast the synchrotron lifetime for Radio particles is 
longer than the age of the nebula, so they fill the entire 
volume. This increase in size at smaller frequencies is 
observed in the Crab Nebula (s^lTllQ]. If one considers 
the pressure anisotropy due to the compressed nebular 
toroidal magnetic field 1^81,(85^, it is also possible to re- 
cover the elongated axisymmetric shape of many PWNe 
(i.e. Crab Nebula, 3C58). 

By comparing observations with the predictions of the 
simple spherically symmetric MHD model it is possible 
to constraint some of the properties of the pulsar wind, at 
least at the distance of the TS. To explain the dynamics 
of the plasma, as well as the emission properties of the 
nebula, the Lorentz factor of the wind is estimated to be 
~ 10^, and the ratio between Poynting flux and kinetic 
energy a ~ 0.003. This shows that nebular properties can 
be used to derive informations on the conditions of the 
pulsar wind at large distances. 

EVOLUTION OF PWNE 

In the analytic model developed by KC84 the SNR has 
only a passive role, providing the confinement of the 
PWN. It is however known that the details of the flow 
structure in the PWN and its emission properties depend 
critically on the boundary condition with the SNR. Given 
the complexity of the interaction and the existences of 
different evolutionary phases, a detailed study of the evo- 
lution of the system PWN-SNR, has been possible only 
recently, thanks to the improvement in computational re- 




FIGURE 1. Left picture: schematic representation of the 
global structure of the PWN in the first phase of its evolu- 
tion inside a SNR. From the center the various regions are: 
1- the relativistic pulsar wind, 2- the hot magnetized bubble 
responsible for the non thermal emission, 3- the free expanding 
ejecta of the SNR, 4- the ISM. These regions are separated by 
discontinuities: a- the wind termination shock, b- the contact 
discontinuity between the hot shocked pulsar material and the 
swept-up SNR ejecta, c-the front shock of the thin shell ex- 
panding into the ejecta, d- the reverse shock of the SNR, e- 
the contact discontinuity separating the ejecta material from 
the compressed ISM, f- the forward SNR shock. Right picture: 
evolution of the PWN size, from free-expansion to sedov phase 
(from Bucciantini et al. 1211 ). 

sources IMHIlillsi. 

It is immediately evident by comparing the energy in 
the SNR (~ 10^^ ergs) to the total energy injected by 
the pulsar during its lifetime (~ 10^^ ergs) that while the 
PWN cannot significantly affect the SNR, the evolution 
of the SNR can have important consequences for the 
PWN. 

A simple ID model for PNWe evolution shows the ex- 
istence of three main phases (for a more complete discus- 
sion of PWN- SNR evolution see Reynolds & Chevalier 
iItsIi and Gaensler & Slane [37]). At the beginning the 
PWN expands inside the cold SN ejecta. The SN ejecta 
are in free expansion, so this phase is somehow referred 
Sisfree expansion phase. During this phase, that lasts for 
about 1000-3000 yr, the pulsar luminosity is high and 
almost constant. This is the present phase of the Crab 
Nebula, and PWNe in this phase are expected to shine 
in high energy X-rays emission. The expansion veloc- 
ity of PWNe in this early stage is known to be a few 
thousands kilometer per second, much higher than typ- 
ical pulsar velocities in the range 50-300 km/s. For this 
reason one can neglect the pulsar kick in modeling young 
objects, and assume the pulsar to be centrally located. As 
the system expands inside the high density, cold, super- 
sonic ejecta of the SNR, a thin shell of compressed ma- 
terial is formed. The evolution of this shell can be easily 
described in the thin-shell approximation, and it is possi- 
ble to derive an analytic self- similar solution L23il . Both 
analytic and numerical results show that the thin shell ac- 
celerates with time. Given that the density of the shell is 
much higher than the enthalpy of the relativistic plasma, 
the shell is subject to Rayleigh-Taylor instability. This is 



supposed to be at the origin of the filamentary network 
of the Crab Nebula H S El • 

During the first phase the PWN is inside the SNR 
shell. The expanding PWN will eventually reach the re- 
verse shock in the SNR shell, which is supposed to re- 
cede to the center of the SNR in a time of order of 
5000-10000 yr. From this moment on the evolution of 
the PWN is modified by the more massive and ener- 
getic SNR shell, and experiences a compression phase 
generally referred as reverberation phase. In the sim- 
ple ID scenario the compressed PWN will then push 
back the ejecta, and the nebula might undergo several 
compressions and rarefactions. However more appropri- 
ate multidimensional studies have shown that the SNR- 
PWN interface is highly Rayleigh-Taylor unstable dur- 
ing compression IIT3I1 . which can cause efficient mixing 
of the pulsar wind material with the SNR. This rever- 
beration phase is supposed to last about 10^ — 10^ years. 
Even if energy injection from the pulsar at these later 
times is negligible, PWNe can still be observed, due to 
the re-energization during compressions. The interaction 
with the reverse shock can lead to a variety of different 
morphological structures if one consider also the pulsar 
proper motion Jsvl Issl Given that the evolution is 
now driven by the SNR reverse shock, the nebula can 
be displaced with respect to the location of the pulsar. 
At the beginning this might result into a system where 
the pulsar is not located at the center of the non ther- 
mal emission, analogous to what is observed in Vela. As 
the system evolves the reverse shock will completely dis- 
place the body of the PWN, creating a relic nebula. The 
relic PWN will mostly contain low energy particles, and 
will be visible in radio, while high energy particles, ob- 
servable in X-rays will be seen only close to the pulsar. 
G24.5 shows indeed this kind of morphology II37I] . At 
later time the SNR starts cooling and eventually the pul- 
sar will become supersonic with respect to the ejecta. In- 
terestingly the location with respect to the SNR where 
this happen does not depend on the pulsar proper motion 
and turns out to be ~ 70% of the radius of the forward 
shock. The pulsar will form around itself a bow- shock 
PWN, and one expects an emission tail to form connect- 
ing the pulsar to the relic PWN. This model apply to the 
morphology and structure of W44 [37]. 

The ultimate phase of a PWN evolution depends on 
the pulsar kick. For slow moving pulsars the PWN will 
expands adiabatically inside the heated SNR, now in 
Sedov phase. Given the absence of energy injection, a 
PWN in this stage is probably only observable as a faint 
extended radio source. Fast moving pulsars can escape 
from the SNR. They will form bow- shock nebulae due 
to the interaction with the ISM, through which they are 
moving at supersonic speeds |[l9l These objects 
are observed both in Ha emission, due to ionization of 
ISM neutral hydrogen, and as long extended cometary- 




FIGURE 2. Upper row: images of PWNe in various evolutionary phases (from Gaensler & Slane Isvh . From Left to right: X-ray 
image of the composite remnant G2 1.5-0.9, free expansion phase; radio image of Vela SNR, displacement of the nebula due to the 
compression of the reverse shock during reverberation; SNR G 327.1-1.1 in radio (red) and X-rays (blue), relic PWN phase; W44 in 
radio, transition to the internal bow-shock phase; PSR B 1957-1-20 in H« (green) an X-rays (red), ISM bow-shock nebula. Numerical 
simulations of the various phases of the PWN-SNR interaction in the case of moving pulsars (from 1,88.. .20.1). Each figure of the 
second row corresponds to systems shown in the first one. 

like source of non thermal emission, due to the shocked 
pulsar wind, now forced to flow in the direction opposite 
to the pulsar motion. 

INNER STRUCTURE 

Recent optical and X-rays images from HST, CHAN- 
DRA and XMM-Newton have triggered new interest of 
the scientific community for PWNe. The new data show 
that the inner region of PWNe is characterized by a com- 
plex axisymmetric structure, generally referred as jet- 
torus structure (Fig. O. This was first observed in Crab 
iisi l92h . and has subsequently been detected in many 
other PWNe S S 112 Izi [3 [fiS S 0^ It 
is characterized by a main emission torus, corresponding 
to what is thought to be the equatorial plane of the pulsar 
rotation. Multiple arcs or rings are often present together 
with a central knot, located in the vicinity of the pulsar, 
and one or two opposite jets along the polar axis, which 
seems to originate close to the pulsar itself. Such struc- 
ture cannot be explained in the simple ID radial model 
by KC84. Even if a higher equatorial energy injection 
[ 16] can qualitatively account for the existence of a main 
torus it is not possible to reproduce the observed lumi- 
nosity using the KC84 model. Shibata et al. |80] was the 
first to point that, the difference in brightness between the 
front and back sides of the torus in Crab Nebula, requires 
a post-shock flow velocity ~ 0.3 — 0.4c, much higher 



than what expected for subsonic expanding flows. More- 
over the existence of an inner ring, separated from the 
torus, contradicted the assumption of a smooth flow from 
the TS. Even more puzzling was the knot which seemed 
to be located inside the wind region. However, the most 
interesting feature was th e je t 16411. b ecause theoretical 
11,^10] and numerical lli [11 SinHl [M] studies of rel- 
ativistic winds from pulsars have shown no presence of 
collimated energetic outflow. Finally photon index maps 
of Crab Nebula |69], and now of Vela [73], show that the 
spectrum flattens moving from the inner ring toward the 
main torus, when a steepening due to synchrotron losses 
is expected. Interestingly the symmetry axis of the jet- 
torus appears to correspond to the major axis of the neb- 
ula, suggesting that the toroidal magnetic field is a key 
element in shaping the inner flow. 

Only recently it was recognized the importance that 
magnetization and energy distribution in the pulsar wind 
have in determining the nebular flow. In contrast with the 
isotropic energy flux assumed in KC84, it was known 
for a long time that the asymptotic solution by Michel 
|[67ll. re cently confirmed with numerical simulations |[l4l 
56, 26], predicts a higher equatorial flux, which naturally 
produces an oblate TS with a cusp in the polar region 
lE[i3]. Given that the obliquity of the TS at higher lat- 
itudes, the post shock flow in the nebula can have speeds 
~ 0.3 — 0.5c. Moreover if hoop-stresses are more effi- 
cient in the mildly relativistic flow, and the collimation 
of a jet might occur in the post shock region I61i,i51j . The 



evident complexity of this scenario clearly prevented any 
sophisticated theoretical model. Only recently, thanks to 
more efficient and robust numerical schemes for rela- 
tivistic MHD [55, ^E^. a detailed numerical descrip- 
tion has been possible. Various numerical results show 
that, based on the theoretical assumption of axisymmet- 
ric energy flux and post-shock hoop stresses, it is pos- 
sible to reproduce the observed features. This picture 
has been subsequently refined by investigating if and 
how the emission features could be used to understand 
the conditions of the pre- shock wind, trying to constrain 
many of the nebular emission properties, like polariza- 
tion and photon index. More recently numerical studies 
have shown that also the short time variability might be 
recovered. 

If one assume for the wind this solution of the force 
free model, the energy flux in the wind has a strong lati- 
tudinal dependence of the form L( ) = ( 1 + a sin ( ) ) , 
where a is a measure of the pole-equator anisotropy, 
while the magnetic field in the wind B{6) oc sm{0). Vari- 
ous numerical simulations of the interaction of such wind 
with the SNR ejecta have been presented |57il33ill8[|341. 
Interestingly these results show that the post shock flow 
is independent on the specific values of Lorentz factor 
or density distribution. In this sense flow dynamics can- 
not be used to constrain the value of the wind Lorentz 
factor or the multiplicity in the wind. In Fig. [3] we can 
see the oblate shape of the TS. Due to the TS shape, the 
flow is slowed down to speed ^ c/3 close to the equa- 
tor, but at higher latitudes, where the shock is oblique, 
the post shock flow is still supersonic. The result of the 
anisotropic energy distribution in the wind is that almost 
all of the downstream plasma is funneled along the equa- 
torial plane. The inner part moves with speed ^ c/3, 
while the outer channels has velocity ~ 0.5c, the value 
expected in order to justify the luminosity distribution in 
the torus of the Crab Nebula LfiO]. 

As the flow expands in this equatorial sheet away from 
the TS, the magnetization increases, until equipartition 
is reached. Given the magnetic field distribution in the 
wind, the flow reaches equipartition close to the equa- 
tor before that at higher latitudes. The post- shock flow is 
only mildly relativistic, so it can be efficiently collimated 
by hoop stresses. Once equipartition is reached, the mag- 
netic pressure prevents further compression of the mag- 
netic field, and the flow is diverted back toward the axis. 
This is the process that causes collimation and the for- 
mation of a jet along the axis itself as seen in Fig. 3. The 
formation of the jet strongly depends on the values of the 
wind magnetization: if it is too small, then equipartition 
is not reached inside the nebula, the equatorial channel 
survives to the edge of the PWN and no jet is formed. At 
higher magnetization equipartition is reached in the close 
vicinity of the TS, and most of the plasma is diverted and 
collimated in a jet. The size and flow velocity in the jet 



are a function of O". For values O" ~ 0.03 the plasma speed 
in the jet is ~ 0.7c in agreement with observation of the 
jet in Crab Nebula f9lS5lE2l]. There is a global merid- 
ional circulation inside the nebula associated with the jet 
formation , with typical speeds ~ 0.1c that might lead to 
local shear instabilities and on the possible mixing with 
cold ions |[66l] . 

It is also possible in numerical models to take into 
account modification of the wind properties associated 
with oblique rotators. While the ener gy d istribution in 
the wind is identical to the aligned case [15, 82], the pres- 
ence of a folded current sheet that might extend to higher 
latitudes, has important consequences on the magnetiza- 
tion. If this striped wind region is dissipated, and this can 
happen either in the wind |65, 54] or at the termination 
shock itself 162, IH], then the magnetization will vanish 
toward the equator and will reach a maximum values at 
intermediate latitudes, depending on the obliquity of the 
pulsar. This unmagnetized region close to the equator, 
adds complexity to the picture presented above. There is 
now an equatorial flow where equipartition is not reaches 
inside the nebula, that survives to much larger distances 
from the termination shock. However equipartition can 
be reached at higher latitudes and the flow is then di- 
verted to the axis and collimated into a jet. Fig. 3 show 
the flow structure. Emission maps based on the results of 
numerical simulations suggest that a striped wind model 
is more promising in explaining the observes inner-ring 
outer-torus structure of many PWNe, while models with- 
out a neutral equatorial region seem more promising in 
the case of single ring nebulae. This show how emission 
features can be used as a probe for the properties of the 
wind. In particular recent results |34] have shown that, 
given the same average magnetization, a bigger striped 
wind region, leads to weaker jets, less intense hoop- 
stresses, and a bigger size of the termination shock. 

Relativistic MHD simulations clearly show the forma- 
tion of a jet, however to properly asses their validity and 
for a detailed comparison with observations it is neces- 
sary to compute the nebular emission. Despite several 
theoretical issues like the particles energy distribution at 
the injection and the efficiency of acceleration, that can 
modify the results, as a first approximation one might 
assume a uniform power-law injection distribution. Un- 
der this assumption, to obtain X-ray maps, one must then 
consider also the adiabatic and synchrotron losses, as 
fluid particles move in the nebula. In Del Zanna et al. ll34ll 
the maximum energy £jnax of the particles' distribution is 
evolved with the flow, accounting both for synchrotron 
and adiabatic losses. In this case the particle distribution 
remains a power low with a break at Emax- A more so- 
phisticated approach would require solving for the entire 
particles distribution at different energies. 

In Fig. 3, we show a CHANDRA image of the Crab 
Nebula, and a map based on a simulation with striped 



FIGURE 3. Fro Left to right: X-ray image of the crab nebula, jet-torus structure. Numerical result of the internal dynamics in 
a PWN, note the oblate shape of the TS, which focus the downstream flow toward the equator, ft33i1 . The flow is then diverted 
back toward the axis but the magnetic hoop stresses, and is collimated into a jet (in this specific case there is an unmagnetized 
region around the equator, and an equatorial floe survives). Simulated X-rays synchrotron map based on numerical simulations and 
simulated map of the photon index from 1. 341 . 



wind in the X-rays. The knot and the inner ring are 
both explained as due to the high velocity flow in the 
immediate post shock region, at intermediate latitudes. 
The main torus observed in the Crab Nebula is due to 
the emission coming from the equatorial channel. Recent 
investigations |91, 34] have shown that the size and shape 
of the striped wind region, have important observational 
consequences. The region corresponding to the inner 
ring might split in several minor arches; the extent of 
the X-rays nebula changes as well as the jet thickness. 
For the case of Crab Nebula simulations also seem to 
disfavor configuration with a higher magnetization close 
to the pole |3]. While modeling of the equatorial region 
does not need any accurate treatment of the synchrotron 
losses (the feature are observed also in radio |[l2n ). these 
are fundamental in understanding the jet, and simulated 
map clearly show that the jet is not visible over the 
background, in optical or radio. 

Interestingly given that the particle energy distribu- 
tion is somehow followed along a streamline, it is pos- 
sible to reconstruct also spectral maps. In Fig. 3 a X-ray 
map of the spectral index based on simulations is pre- 
sented. It reproduces the main observed properties of the 
Crab Nebula |69]. We also see that simulations produce 
maps where the spectrum appears to flatten moving to- 
ward the main torus, without the need to assume any 
re-acceleration, this agree also with recent results about 
Vela 1 73]. The reason is the Doppler boosting effect: for 
low speeds, observing at a given frequency implies sam- 
pling particle all at the same energy (if one neglects vari- 
ations of magnetic field). When velocities are relativistic, 
the energy of particle responsible for emission at a given 
frequency, depends also on their Doppler boosting. The 
flatter spectrum in the torus, is due to the higher speed. 
There are still problems to recover the correct spectrum 
in the jet: the simulated spectrum is much steeper than 



what is observed, due to excessive synchrotron losses. To 
reproduce observations, one must assume some form of 
dissipation and re-energization along the axis. This could 
be associated with local instabilities of the toroidal mag- 
netic field |i6J which are not captured by axisymmetric 
simulations, but for which there are many observational 
evidences Jv^, IzS Esl EH] • Integrated spectra can also be 
used to constrain the average nebular magnetization. Pre- 
liminary results show that O" > 0.03 is required and con- 
figuration with a smaller striped wind region are favored. 
Interestingly integrated spectra show a flattening above 
the CHANDRA band, which has been interpreted as due 
to the high velocity in the posts-hock region. This could 
be a numerical artifact, or it may hint to the possibility of 
an injection break at high energies (either associated with 
acceleration process itself or with a changing efficiency 
along different portion of the TS). Further investigation 
will be necessary to understand this effect. 

Polarization might also be used to derive useful infor- 
mation about the magnetic field structure. While emis- 
sion maps are strongly dependent of the flow velocity 
inside the nebula, they are not much sensitive to the 
presence of a small scale disordered component of the 
magnetic field. In this case polarization might be use- 
ful. Given that only optical polarization is available, any 
study should be limited to the brightest features, like the 
wisps, inner ring and torus. The effect of flow velocity 
on the polarization angle has been discussed by Buc- 
ciantini et al. |25] and Del Zanna et al. ll34|] By com- 
parison with old optical polarization map 179*], one can 
immediately see that some of the properties are recov- 
ered, among whom a region which can be identified with 
the bended jet of Crab. More recent data (Graham, pri- 
vate communication) have also shown that the knot has 
a high degree of polarization, and that the polarization 
angle is consistent with the nebular origin of this feature. 



FIGURE 4. Variability in the wips region due to ion compression | 83]. The pictures represent the X-ray luminosity at three 
different instants. Notice the bright features moving outward. 
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FIGURE 5. The Plot show the X-ray luminosity, based on 
MHD simulations, along the axis of a PWN as a function of 
time. Note the 2-year quasi periodicity in the wisps region and 
for the knot. Features appear to move out from the inner TS 
region, the propagation velocity is ~ 0.5c in the interior and 
slows down at larger distances. 

TIME VARIABILITY 

The former discussion has focused on the explanation 
of the main observed features, that appear to be quite 
persistent on long time- scale. However it is known that 
PWNe show a short time variability at high energies. 
Variability of the wisps in the Crab Nebula has been 
known since the first high resolution observations llisi 
fl2n The jet in Vela appears to be strongly variables 
and there is evidence for variability of the main 



rings ll73n . Variability is also observed in the jet of Crab 
||m[59[r and B 1509 Isi]]. B 1509 shows also variability 
in the inner ring. Variability in the jet, which usually 
has time-scale of years, has been associated to kink or 
sausage modes in the strong toroidal filed, or even to 
fire-hose instability 1 84] in the case of Vela. On the other 
hand the wisps show variability on time-scale of months, 
much shorter than typical sound crossing time, in the 
form of an outgoing wave pattern, and recent results for 
the Crab Nebula suggest the presence of a year-long duty 
cycle. 

moreFor a long time the only model capable of repro- 
ducing the observed variability was the one proposed by 
Spitkovsky & Arons |83]. The model is based on the as- 
sumption that ions (or high energy electrons) are present 
in the equatorial region of the wind. The presence of ions 
is consistent with the fact that, in kinetic simulations of 
acceleration in a strong shock, a pure pair plasma do not 
produce any power-law distribution |2]. Given that ions 
have a much larger larmor radius, of order of the size 
of the wisp region, they introduce a substantial deviation 
from a pure fluid picture: electrons are compresses by 
the gyration of ions and emission is enhanced. By select- 
ing the fraction of energy in the ion component one can 
reproduce the observed time scale variability and the av- 
erage distance between the torus and inner wisps. How- 
ever results are based on a simplified ID model that does 
not take into account the energy flux anisotropy, and the 
shape of the TS. 

MHD nebular models however show the presence of 
high velocity flow channels inside the nebula, and it is 
possible that shear instabilities can be at the origin of 
the observed variability |27l |7t]. It was already noted in 
simulations by Bogovalov et al. fl^] that the synchrotron 
emissivity inside the nebula varies. Results by Komis- 
sarov & Lyubarsky [57] also suggest that the flow might 
have a feedback action on the TS, causing it to change 



shape, and thus inducing a change in the wisps, which 
originate from posts-hock flow. It was not clear however 
if the observed variabihty in the internal flow pattern, 
could result in the observed variability of the wisp, if 
a wave pattern was recovered, and what were the time- 
scales. Recent results show that MHD models can ac- 
count for the observed variability. The shear between the 
outgoing equatorial channel and the back-flow that is di- 
verted by hoop stresses toward the axis, causes the for- 
mation of eddies on scales typical of the TS, that are 
subsequently advected away from the TS. These eddies 
in turn change the shape of the TS, and the Doppler 
boosting responsible for the arches and rings. Simula- 
tions show that there is a typical duty-cycle of about 1 
year, and that the outgoing wave pattern is recovered. 
Typical speeds for such waves features is also in agree- 
ment with observations. 



GAMMA RAYS 

Recently the interest of the scientific community has 
focused on the new gamma-ray results L42, 76]. HESS 
has been able in a few cases to resolve some PWNe both 
young ones in the free expansion phase and older one 
undergoing reverberation [ 30, 381 • New results will soon 
come with GLAST which will be able to investigate also 
the high energy MeV synchrotron emission from the tail 
of the particle energy distribution. Even if this results 
lack the high spatial resolution of the X-ray images they 
can provide useful information about PWNe. 

The emission at MeV energy observed from young 
object will enable us to put constraints on the acceler- 
ation mechanism. Interestingly, given the variability of 
the wisps structure, one might expect that the emission 
at high energy should show similar variability, on compa- 
rable time- scales. On the other hand given the short life- 
time for synchrotron losses of these particles, they might 
can be used to investigate short time variability in the 
wind (X-ray particles have longer lifetime). In principle, 
for example in the case of nulling pulsar [58], one could 
use the emission in this band to understand how the wind 
is modified in the null phase. 

The emission at GeV energies is assumed to be from 
Inverse Compton scattering on background radiation, 
and in the case of young objects also the self synchrotron. 
Interestingly the spectral properties of the comptonized 
radiation can be used to derive information about parti- 
cles which are supposed to emit synchrotron in the UV. 
Standard models for the evolution of the energy distribu- 
tion function inside PWNe predict that the synchrotron 
break in the distribution should be in the region corre- 
sponding to these UV emitting particles. GeV emission 
could in principle help provide an independent constraint 
on the magnetization in the nebula. Moreover if protons 



are present in the pulsar wind, then it is possible that they 
might contribute to the gamma ray emission vis p-p scat- 
tering and pion decay fT, 48]. 

Presently, evolutionary models for the high energy 
emission from PWNe, have been developed only in the 
case of young system during the free expansion. How- 
ever, many of the PWN observed in the GeV band, are 
old objects in the reverberation phase, and it is important 
for a correct modeling of the emission to take into ac- 
count the effect of the compression. This can lead to sig- 
nificant modifications both in the case of a pure leptonic 
model, due to the re-compression of the nebula, and also 
in the case of adronic emission, given that an efficient 
mixing with the colder SNR material is expected (thus 
increasing the density of target protons). 

In many of these objects the bulk of the gamma emis- 
sion is not centered on the pulsar, and the displacement 
is too large to be explained in term of a moving pulsar 
leaving behind a relic PWN. Two possible explanations 
have been invoked: either a SNR exploding in a density 
gradient |TH], or a bow-shock. In the first case the re- 
verse shock will collapse off-center with respect to the 
SNR, displacing the PWN. Observationally one should 
look for evidences of density gradients i he ISM, for ex- 
ample in the form of molecular clouds or star forming 
regions. In the second case, simulations and observations 
suggest that the pulsar wind back flow in the bow- shock 
tail can have high speeds up to fraction of the speed of 
light|24, 50]. Particle responsible for the inverse comp- 
ton emission in the gamma-ray band can be advected to 
large distances from the pulsar. In this case the displace- 
ment of the emission with respect to the pulsar should be 
aligned with the pulsar proper motion. 

CONCLUSION 

In the last few years, the combination of high resolution 
observations, and numerical simulations, has improved 
our understanding of the evolution and internal dynamics 
of PWNe. We can reproduce the observed jet- torus struc- 
ture and we can relate the formation of the jet in the post 
shock flow to the wind magnetization. Simulated maps 
can reproduce many of the observed features, including 
the details of spectral properties. Results suggest that, the 
best agreement is achieved in the case of a striped wind, 
even if MHD simulations are not able to distinguish be- 
tween dissipation of the current sheet in the wind or at the 
TS. Results also suggest that it is possible to use X-rays 
imaging to constrain the pulsar wind properties; already 
the rings and tori observed in many PWNe have been 
used to determine the spin axis of the pulsar Inter- 
estingly in Crab the inner ring does not appear boosted, 
while the wisps (which are interpreted as its optical coun- 
terpart) are. 



There are still however unsolved questions, and pos- 
sible future developments for research in this field. All 
present simulations are axisymmetric, and none is able to 
address the problem of the stability of the toroidal field, 
and cannot reproduce the observed variability in the jet. 
It is not clear if small scale disordered field is present in 
the inner region (may be a residual of the dissipation in 
the TS of the striped wind). A combination of simula- 
tions and polarimetry might help to answer this question. 
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